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Abstract 

In this paper we first obtain the order of stress singularity for a dynamically propagating self-affine fractal 
crack. We then show that there is always an upper bound to roughness, i.e. a propagating fractal crack 
reaches a terminal roughness. We then study the phenomenon of reaching a terminal velocity. Assuming 
that propagation of a fractal crack is discrete, we predict its terminal velocity using an asymptotic energy 
balance argument. In particular, we show that the limiting crack speed is a material-dependent fraction of 
the corresponding Rayleigh wave speed. 
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1 Introduction 

A theoretical framework for including inertial effects during a rapid crack growth was first proposed by Mott 
[1948], who adopted the analysis of Griffith [1921] as a starting point. The idea is based on a simple addition of a 
kinetic energy term to the expression for the total energy of the cracked system. According to Mott's extension of 
Griffith's criterion, the requirement that the system remains in thermodynamic equilibrium with its surroundings 
as the crack extends leads to the following expression in terms of the well-known fracture parameters: G — 27 = 
dT/da, where G is energy release rate, 7 is specific surface energy T is kinetic energy density, and "a" is the 
characteristic length of the crack. Mott defined a domain R that receives stress- wave "messages" from the crack 
tip and then argued that the total kinetic energy can be written asT= \pv 2 J R [(du x /da) 2 + (du y /da) 2 ]dxdy. 
While Mott's analysis lacks rigour, it is instructive in the way it highlights some of the important features of a 
running crack without excessive mathematical complication. 

The first important contribution to the problem of a moving crack with constant velocity was the work of 
Yoffc [1951]. The Yoffc problem consists of a mode I crack of fixed length traveling through an elastic body at a 
constant speed under the action of uniform remote tensile loading. Yoffc [1951] obtained the stress distribution 
near the tip of a rapidly propagating crack in a plate of isotropic elastic medium. The result was that the 
stresses depend on the crack tip velocity and reduce to the solution of Inglis [1913] when the velocity is zero. 

Roberts and Wells [1954] used Mott's extension of Griffith's criterion to predict the limiting velocity of the 
crack extension. By taking the boundary of the region R to be a circle of radius r centered at the crack tip, 
they estimated r ps cot, where Co = \J E/ p is the longitudinal sound wave speed. They defined this cutoff region 
as the border of the disturbed zone by the stress waves emanated from the crack tip. Using this assumption 
and taking a stress field similar to that of the static case they roughly estimated the limiting crack velocity to 
be about 0.38c when v = 0.25. Steverding and Lehnigk [1970] studied the problem of the response of cracks to 
stress pulses and found an equation of motion for such cracks. They also obtained the limiting velocity of crack 
extension caused by stress pulses by using asymptotic solutions to be about 0.52cr, where cr is the Rayleigh wave 
speed. There have also been some other efforts on finding the equation of motion for dynamically propagating 
cracks. Berry [1960a,b] was the first to find an equation of motion for dynamic propagation of cracks. He found 
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out that the details of the motion of a crack are determined by the state of stress at the point of fracture, 
and that the observed critical stress is (infinitesimally) greater than that given by the Griffith's criterion and is 
probably determined by the size of the defect in the sample and the rate of straining. He also obtained solutions 
for the equation of motion for fracture in tension and fracture in cleavage in both constant force and constant 
velocity cases. See Bouchbinder, et al. [2010] for a recent review of dynamic fracture mechanics. 

The inadequacy of the classical fracture mechanics theories in problems such as predicting infinite strength 
for elastic bodies without any cracks, for example, was the motivation for some researchers to propose new 
failure theories. Novozhilov [1969] introduced a non-local stress criterion and gave the condition of the brittle 
crack propagation in mode I as a* = {<J v {x))q = a c , where (.)g° is spatial averaging over the interval [0, do], 
cr y (x) is the complete (not only asymptotic) stress field around the crack tip (x = 0), er c is the ideal strength of 
the material, and ao is the minimum admissible crack advance named by him a fracture quantum. According 
to Novozhilov this criterion can be used only with the complete expression of the stress-field, and not with 
its asymptotic form. However, the complete expression is rarely known. Another restriction in Novozhilov's 
approach was that the size of fracture quantum assumed to be the atomic spacing. Pugno and Ruoff [2004] 
introduced their so-called quantized fracture mechanics (QFM) approach, which modified Novozhilov's theory. 
In QFM, the restrictions of Novozhilov's theory were removed and this made QFM a useful approach for analysis 
of very short cracks (see also Krasovs'kyi [2006]; Morozov and Petrov [2002]; Cornetti, et al. [2006]; Lcguillon 
[2002] for more related works). In their approach the differentials in Griffith's criterion were replaced by finite 
differences (see Wnuk and Yavari [2008] for a discussion). For vanishing crack length, QFM predicts a finite 
ideal strength in agreement with the prediction of Orowan [1955]. 

In most models in fracture mechanics cracks are assumed to be smooth for mathematical convenience. 
However, in reality fracture surfaces are rough and "roughness" evolves in the process of crack propagation. 
Fracture surfaces of many materials of interest are fractals, a fact that has been experimentally established 
by many researchers. A fractal dimension (roughness exponent) is not enough to uniquely specify a fractal 
set and this is why all one can hope for achieving having only a fractal dimension (roughness exponent) is 
a qualitative analysis. Effects of fractality on fracture characteristics of rough cracks have been investigated 
by several groups in the past two decades (see Mosolov [1991]; Gol'dshtein and Mosolov [1991, 1992]; Balankin 
[1997]; Borodich [1997]; Carpinteri [1994], Cherepanov, et al. [1995]; Xie [1989, 1995]; Yavari, et al. [2000]; Yavari 
[2002]; Yavari, et al. [2002a,b]; Wnuk and Yavari [2003, 2005, 2008, 2009], and references therein). Here our 
interest is to estimate the observed terminal velocity of a rough crack propagating dynamically in an clastic 
medium. 

Wnuk and Yavari [2008] extended quantized (finite) fracture mechanics ideas for fractal cracks. They pre- 
sented a modification of the classical theory of brittle fracture of solids by relating discrete nature of crack 
propagation to the fractal geometry of the crack. Their work is based on the idea of using an equivalent smooth 
blunt crack with a finite radius of curvature at its tips for a given fractal crack [Wnuk and Yavari, 2003, 2005]. 
By taking the radius of curvature of the equivalent blunt crack as a material property, they showed that fractal 
dimension of the crack trajectory is a monotonically increasing function of the nominal crack length. This result 
was an analytical demonstration of the mirror-mist-hackle phenomenon for rough cracks. Later they showed 
that assuming a cohesive zone ahead of a fractal crack, the size of the cohesive zone increases while the crack 
propagates [Wnuk and Yavari, 2009]. 

To our best knowledge, the only contributions related to the dynamic fracture of fractal cracks are Xie [1995] 
and Alvcs [2005]. Xie [1995] introduced a fractal kinking model of the crack extension path to describe irregular 
crack growth. Then by using the formula proposed by Freund [1998] for calculating dynamic stress intensity 
factor for arbitrary crack tip motion he calculated the stress intensity factor for the assumed fractal crack path. 
He concluded that the reason for having terminal velocities lower than the Rayleigh wave speed is the fractality 
of the crack path. Alves [2005] used a fractal model for rough crack surfaces in brittle materials. He tried to 
explain the effect of fractality of fracture surfaces on the stable (quasi-static) and unstable (dynamic) fracture 
resistance. He concluded that fractal dimension has a strong influence on the rising of the R-curve in brittle 
materials. He also argued that the reason for having terminal velocities lower than Rayleigh wave speed is the 
roughness of the fracture surfaces that makes the nominal (projected) and local crack tip velocities different. 

Dulancy and Brace [1960] modified Mott's analysis of energy balance and showed that for a crack of initial 
and current lengths cq and c, respectively, c = «l(1 — co/c) compared to the similar equation c = vl(1 — cq/c)? 
obtained by Mott. Here vl is the terminal velocity and is proportional to y/E/p. They also carried out tests of 
terminal velocity on PMMA specimens and showed that the measured velocity differs only about ten percent 
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from the predicted value by Roberts and Wells [1954]. Recently, Chekunaev and Kaplan [2008, 2009] studied 
the terminal velocity by replacing a mode I crack under pressure on its faces by considering cohesive forces 
and replacing the crack by an equivalent distribution of dislocations. They obtained simple expressions for the 
potential and kinetic energies of the environment of the moving crack. They also obtained an expression for 
an equivalent mass for the crack tip, i.e. a point mass that has the same kinetic energy as the whole cracked 
system. Their equivalent mass depends on a truncation radius R in the form of ln(2i?/a), where a is the half 
crack length. For a uniform external pressure po they showed that the crack tip speed can be expressed as 
i> = vl(1 — a cr /a), where a cr depends on both mechanical properties and pq. Their terminal velocity has the 
form Vl = g(v, Rj a)\J E j p that they approximate by Uj, = giv)\J ' E / p, for functions g and g that are given 
in [Chekunaev and Kaplan, 2008, 2009]. It should be mentioned that their equivalent mass is positive for all 
R/a > and v e (—1.0,0.5). Note that their equivalent mass is independent of the crack tip speed and hence 
kinetic energy is an increasing function of crack tip velocity. 

This paper is organized as follows. In §2 an asymptotic method will be used to determine the order of 
stress singularity for a dynamically propagating fractal crack. In §3 dynamic propagation of a fractal crack is 
investigated. We first show that in the intermediate crack growth regime, i.e. after the initial phase of crack 
growth, roughness change is so small that a terminal roughness exponent can be assumed. The phenomenon of 
reaching a limiting speed is predicted using some simplifying assumptions. The predicted limiting crack speeds 
for different brittle amorphous materials are shown to have good agreement with the experimental results. 
Finally, conclusions are given in §4. 



2 Order of stress singularity for a dynamically propagating fractal 
crack 

Consider a crack, whose tips are growing in opposite directions with equal velocities. 1 The crack unloads some 
area of the body, and while propagating, size of the unloaded area will increase. The form of the stress field in 
the close vicinity of the crack tip is of interest as almost all fracture energy will be consumed through different 
processes in this zone. In the following, the asymptotic behavior of some important parameters that contribute 
to energy balance will be investigated. Here, we assume that the singularity of stress is of the form r - ' 3 , where 
r is distance from the crack tip. In the sequel we find an expression for the order of stress singularity 

Several experimental observations confirm that energy consuming phenomena such as temperature rise, 
acoustic and phonon emissions, etc. occur in the close vicinity of the moving crack tip. There is one common 
aspect in all these phenomena; they all have a kinetic origin. To be more precise they are all results of fast 
movements or oscillations of the particles around the moving crack tip. This means that the velocity of the 
particles around the crack tip is of great importance. The particle velocity at a point depends on two main 
parameters: the crack tip speed and the strain at the point [Freund, 1998]. Asymptotic behavior of particle 
velocity in the close vicinity of a smooth moving crack tip is similar to that of stress for all modes of fracture. 
For the classic case we have iii ~ r~2 . The classic case is a limiting case of fractal model, which is when the 
roughness (Hurst) exponent (H) of a self-affinc crack trajectory is equal to unity. In the case of a fractal crack, 
we know that u; = fi(v, Kj, r, E) and hence dimensional analysis tells us that 

u, = «-^* 4 (M), (2-1) 

where vpj is a dimcnsionless function and v is short for ^nominal- 

Asymptotic behavior of the true length of the crack trajectory is also important. The first experimental 
study on fractal characteristics of fracture surfaces was carried out by Mandelbrot, ct al. [1984] who showed 
that the fracture surfaces of steel are fractals. Since then many experimental investigations have been done. For 
example, the investigations on the concrete fracture surfaces by Saouma, et al. [1990] and Saouma and Barton 
[1994] showed that the fracture surfaces of concrete are also fractals. Based on the experimental observations, 
we assume that fracture surfaces are self-affine fractals. The asymptotic behavior of the true crack length growth 



1 To avoid problems with stress wave reflections we can assume that the crack is semi-infinite. This assumption will not change 
anything in the following analysis. 
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AL is [Mandelbrot, 1985, 1986a,b; Yavari, et al., 2002a] 

AL ~ (2.2) 

where £ is some characteristic crack growth length. 2 . 

The changes of kinetic and strain energies of the body 1Z when the crack length growth is AL arc (summation 
over repeated indices is implied): 



1 , . , . „ fa 



T 2 



AT = I -pu.ii.dA and AU e = / — dA, (2.3) 
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where dA is the area element. If it is assumed that the change of the strain and kinetic energies is dominant 
in a small neighborhood 1Z S around the crack tip, in the above relations 1Z can be replaced by 1Z S . Kinetic and 
strain energy changes have the following asymptotic expressions: 

AT~f- 2f3 and AU e ~ £ 2 ~ 2 ' 3 . (2.4) 

In addition to the kinetic and strain energies, there is another important energy term, namely the surface energy. 
By assuming a constant specific surface energy per unit of a fractal measure, the surface energy that is required 
for the formation of a self-affine fractal crack has the following asymptotic behavior 3 : 

AU S ~£&. (2.5) 

For quasi-static crack growth, Griffith's criterion can be written as AU e + AU S = 0, while for dynamics crack 
growth it is written as AU e + AU S + AT = 0, where AU e is the change of the strain energy in the body due to 
crack growth, AU S is the required energy for the formation of the new fracture surfaces, and AT is the change 
of kinetic energy in the body. Using the equality of asymptotic expressions of energy terms the order of stress 
singularity can be obtained as 

£ 2 -^~ei. (2.6) 

Thus 

which is identical to that of a stationary crack [Mosolov, 1991; Gol'dshtein and Mosolov, 1991, 1992; Balankin, 
1997; Yavari, et al., 2000, 2002a; Yavari, 2002]. Therefore, the stress field has the following asymptotic form: 



r 




where 0=1 2H ' 2 < H < ? . (2.8) 

< H < \ y ' 



3 Terminal velocity of rough crack growth in brittle materials 

From many studies of fracture surfaces formed in brittle materials, it is believed that the surfaces created by the 
process of dynamic fracture have a characteristic structure, referred to as mirror-mist-hackle in the literature. 
This structure has been observed to occur in materials as diverse as glass and ceramics, noncrosslinkcd glassy 
polymers such as PMMA and crosslinked glassy polymers such as Homalite 100, polystyrene and epoxies (for 
more details see Lawn [1993]; Gao [1993]; Hauch and Marder [1998]; Ravi-Chandar [1998]; Fineberg and Marder 
[1999] and references therein). 

The prediction and measurement of the crack tip speed has received great attention from researchers in 
the field of dynamic fracture. As was mentioned earlier, Roberts and Wells [1954] were the first to find a 

2 There are different definitions of fractal dimension, e.g. box dimension (Dg), compass dimension (Bp), and mass dimension 
{Dm)- I n the case of a self-similar fractal, all of these dimensions have the same value, but this is not the case for sclf-affinc fractal 
sets. The local values of the box dimension (using small boxes) and mass dimension (using small radii) are both 2 — H. The 
compass dimension has a local value j^. For more details see [Mandelbrot, 1985, 1986a, b]. 

3 For every short cracks surface energy is length dependent [Ippolito, et al., 2006]. However, we are interested in obtaining the 
limiting crack velocity that corresponds to crack lengths much larger than any fracture quantuum and hence specific surface energy 
is a material constant. 
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theoretical prediction for limiting crack tip velocity. Their calculations based on the Mott's extension of Griffith's 
criterion predicted crack tip speed of 0.38co (for v = 0.25), where c = \J Ej p. Steverding and Lehnigk [1970] 
also predicted the terminal velocity by using an asymptotic solution and found it to be about 0.52cr. Many 
researchers argue that the maximum velocity attainable by any moving surface of discontinuity should be 
identified with the velocity of the Rayleigh surface waves (cr). Ravi-Chandar [1998] reached the following three 
major conclusions about the crack speed measurements: i) There is an upper limit to the speed with which 
dynamic cracks propagate, ii) This limiting crack speed is significantly lower than the Rayleigh surface wave 
speed of the material, iii) The limiting speed is not a fixed fraction of the Rayleigh wave speed; this fraction is 
material dependent. The data gathered by Fineberg and Marder [1999] indicates that in amorphous materials 
such as PMMA and glass, the maximum observed velocity of crack propagation barely exceeds about 1/2 of the 
predicted value. 4 In the following we predict the terminal crack tip velocity using an asymptotic energy balance 
argument. 

To illustrate the process of reaching a constant crack tip speed, suppose that an infinite domain 1Z with 
an initially smooth crack is subjected to remote tensile stresses (Too. The crack unloads some area of the body 
1Z C that can be approximated by a disk 1Z S of radius r s [Yavari, et al., 2002a]. To specify this circle we need 
to define a characteristic length for the problem. There is experimental evidence that the dynamic fracture 
processes approach a steady state and thus taking the crack length as a characteristic length will contradict 
a steady state condition. 5 Therefore, we seek a new characteristic length in the problem. Here the fracture 
quantum is taken as the material characteristic length. Therefore, for a fractal crack the radius of the disk 1Z S 
is assumed to be proportional to ao, i.e. r s ~ ao from dimensional analysis arguments. More precisely, for a 
sclf-affine crack, the following relation holds for the radius of the dominance region of strain energy release: 

r s = a $(H), (3.1) 

where $ is a dimcnsionless function. We assume that the crack is initially smooth (H = 1). At time t = the 
strain energy density reaches a critical value and suddenly the crack starts to grow. As the nominal length of 
the crack increases, the roughness of the fracture surfaces increases as well (H decreases) due to the mirror-mist- 
hacklc transition phenomenon. We postulate that crack surfaces reach a terminal roughness. It should be noted 
that a self-affine fractal model for H < \ is a plane-filling set and hence the limiting roughness lies in the range 
i < Hl, < 1. Let us justify our postulate of reaching a terminal velocity using a crack branching argument. 
The required surface energy for the formation of a self-affine fractal crack has the asymptotic behavior £tt . To 
estimate the actual growth of a self-affine fractal crack, suppose that the nominal crack growth step is equal to 
nao, where ao is fracture quantum and n > 1. According to fractal geometry concepts the actual growth length 
AL of a self-affine fractal crack with the nominal growth step size of nao has the asymptotic form AL ~ aon^ . 
Therefore, the required surface energy has the asymptotic behavior of the form n«. Now we argue that for 
some values of H the required surface energy for the formation of two new surfaces (assuming that these new 
cracks are initially smooth) will become smaller than the required surface energy for the continuation of the 
single (roughened) crack. 6 For roughness exponents smaller than this limiting roughness (denoted by Hl) one 
can write: 

(ti)tt > 2{n) for H < H L - (3.2) 

Therefore, it is probable that by reaching Hl the increase of energy flow toward the crack tip causes branching 
of the crack. This roughness limit can be estimated by solving n~» = 2n for different values of n. 7 The results 
are presented in Fig. 3.1. Note that for each n the acceptable values of H for having a single crack are H > Hl- 
There are different arguments in the literature for explaining the branching phenomena in fracture but many 

4 Note that this is true only for isotropic materials. In an anisotropic body, terminal velocity can reach up to ninety percent 
of the Rayleigh wave speed. See the review by Fineberg and Marder [1999] for more details. Here we restrict ourselves to crack 
propagation in an isotropic medium. 

5 In the case of a semi-infinite crack there is no characteristic length. 

6 Eshelby [1971] suggested that a crack would branch when the energy going into the creation of a single propagating crack is 
enough to support two single cracks. For more details see Fineberg and Marder [1999]. 

7 Note that the equation that we need to solve is k{H)n~H = 2n for a function k such that fc(l) = 1. As we do not have an 
explicit form for k(H), we assume it is approximately equal to unity. If we assume other constant values for k(H) the only change 
will be value of the approximate terminal roughness. It seems that the choice k = 1 leads to a reasonable terminal roughness in 
agreement with experiments. We should also emphasize that the exact value of this terminal roughness will not change any of the 
subsequent results. 
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of them are not in agreement with the experimental results. 8 Note that in the range n £ [10, 100], Hl ~ 0.8 
is almost constant. For large n, H increases but very slowly. Note that as n — > oo, Inn = o(n c ) V e > 
[Blcistein and Handclsman, 1986]. This means that Inn increases indefinitely but very slowly as n increases. 
Note also that for short cracks surface energy depends on n and hence in Fig. 3.1, n < 5 is not shown. 

1.0 

0.9 
H L 0.8 
0.7 
0.6 



10 20 30 40 50 

n 

Figure 3.1: Limiting roughness for different values of n = AL nornina i/aQ. 

Since the pioneering work of Mandelbrot, ct al. [1984] understanding the morphology of fracture surfaces has 
been a very active field of research. Much recent effort has been focused on characterizing fracture surfaces in 
terms of a roughness exponent. There are many researchers who argue that there exists a universal roughness for 
fracture surfaces. Their studies on the fracture surfaces of different materials indicate that for both quasi-static 
and dynamic fracture a universal roughness exponent of approximate value 0.8 can be obtained for values of 
£ (scale of observation) greater than a material-dependent scale, £ c [Bouchaud, et al., 1990; Bouchaud, 1997, 
2003; Mal0y, et al., 1992; Daguier, et al., 1996, 1997; Ponson, et. al, 2006]. What we see from our simple crack 
branching argument is in agreement with these experimental studies. 

By taking the nominal crack growth step to be nao the actual crack growth length for a self-afiine crack 
propagation can be estimated as follows: 

AL actual = \j Sx 2 + 5y 2 ~ a \] n 2 + ni ~ a nV (n > 1). (3.3) 

Note that we assume that the crack lies in the (x, y) plane with nominal growth in the x direction. Let us define 

i 

the nominal crack-tip velocity by f nominal = AL '^° inal = and the actual velocity by factual = " a°° , where 
At is the required time for growth of the fractal crack by the nominal amount of nao- Therefore, the following 
relation holds between nominal and actual velocities of the crack tip: 

factual ~ ni^ 1 V nominal- (3.4) 

l i 

When the crack roughness reaches its limiting value = H^n), n H L ( - n) = 2 and hence factual ~ 2f n0 minai- 
Now if the limit of factual is Cr, we see that limit of f nom inai is about \cr from this rough estimate. Gao [1993] 
used a wavy-crack model and observed that depending on roughness local crack tip velocity can be as large as 
twice the apparent crack tip speed. He also observed that when apparent crack speed is \cr dynamic energy 
release rate is maximized. 

8 Once a crack bifurcates, single crack models are, of course, no longer valid. Therefore, a theory describing a single crack 
can, at best, provide a criterion for when crack branching occurs. A number of such criteria for the onset of crack branching 
have been proposed. The criterion due to Yoffe [1951] and extremal energy density criteria [Sin, 1973; Theocaris and Georgiadis, 
1985; Ramulu and Kobayashi, 1986; Adda-Bcdia and Aniar, 1996] all suffer from the common problem that the velocities pre- 
dicted for the onset of branching are much higher than the observed velocities in the experiments. Additional criteria such as 
postulating a critical value of the stress intensity factor, have not been consistent with experiments [Ramulu and Kobayashi, 1985; 
Arakawa and Takahashi, 1991] since measurements at the point of branching show considerable variation of the stress intensity 
factor Kj. 
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Reaching the terminal roughness exponent Hl has another important consequence; after reaching Hl the 
radius of the dominance zone of strain energy release will remain unchanged, i.e. r s reaches a constant value. 
Now we are back to the main problem of having a terminal velocity for crack propagation in brittle materials. As 
we concluded earlier from the experimental observations, almost all the released strain energy will be converted 
into kinetic energy. Therefore, we are concerned with the changes of strain and kinetic energies, i.e. AU e and 
AT. The dominant change of strain energy can be written as follows: 

AU e = / ^cTijeijdA. (3.5) 

Similar to the case of dynamic fracture of a smooth crack, for a fractal crack the following stress and strain 
fields at the moving crack tip are assumed: 

(Jijir, 9, v) = K{{v)r-^ l3 {e, v) and e. y -(r, 9, v) = K[{v)r^ ] C m Vkl(0, »)• ( 3 - 6 ) 
Substituting the above asymptotic fields into Eq.(3.5), we obtain: 

1 r 



K{{v) r~ 2 ^ tJ (9,v)C ljkl ^ k i(9,v)dA. (3.7) 



In the case of a mode I smooth crack, stress field explicitly depends on the instantaneous crack tip speed v(t) 
[Freund, 1998]. An immediate consequence of this is that the near tip field for the time-dependent motion is 
identical to that for steady state crack growth in the same material up to a time-dependent proportionality 
constant. This was demonstrated by Freund and Clifton [1974], Nilsson [1974], and Achenbach and Bazant 
[1975], each of whom compared asymptotic solutions for nonuniform crack growth with earlier results based on 
the assumption of steady growth obtained by Cotterell [1964], Rice [1968], and Sih [1970]. 

Freund [1972a, b] suggested an indirect method for determining the dynamic stress intensity factor for mode 
I crack propagation with nonuniform crack growth and a general loading condition. The result was that the 
dynamic stress intensity factor for arbitrary motion of the crack tip is proportional to the corresponding quasi- 
static stress intensity factor with a universal proportionality constant, i.e. Ki(a,v) = k(v)Ki(a,0), where a is 
the instantaneous crack characteristic length, and 

«„) = )-"'" . (3.8) 

S +( l / V )\f l - V / C <1 

where cr, Cd, and v are the Rayleigh wave speed, dilatational wave speed, and the crack tip velocity, respectively, 
and S+(l/v) is close to unity. Therefore, k(v) can be approximated by 

k(v) = , (3.9) 

where £ = Cd/cR. 9 We assume that the above result holds for a fractal crack, i.e. 10 

K{{a, v) = k(v)K{(a, 0). (3.10) 

The quasi-static stress intensity factor for a fractal crack of projected length 2a subjected to uniform far-field 
stress a°° is K f IQ = ^(^^Aaf^P [Yavari, 2002; Wnuk and Yavari, 2003]. Thus, for very long cracks 
(intermediate crack growth regime) the rate of change of stress intensity factor becomes vanishingly small, i.e. 

^-—^-0 as a^oo. (3.11) 
oa 

Now, we can use the above arguments in the calculation of the strain energy release. If roughness reaches 
its terminal value Hl-, the size of the dominance zone of strain energy release 1Z S will became approximately 
constant and as a result the order of stress singularity will remain unchanged. In addition to this, increasing 



9 Note that for plane strain c^/cr € [1.635, oo], and for plane stress c^/cr £ [1.635, 2.145]. 
10 Xie [1995] made a similar assumption for each prefractal crack trajectory in his fractal model. 
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the nominal crack length to large values the change in the stress intensity factor becomes vanishingly small. In 
other words, the stress intensity factor becomes approximately constant. Under these conditions, strain energy 
will change only due to the change of the nominal velocity of the crack tip. Strain energy in the disk 1Z S is 
(F->J?i,o»l): 



AU e (v) ~ \k{vf 



I 


-"■JO 


ln B {H L ) 





K 



JO 



' k(v) 2 H L rT L 



r~ (0, v)C m Z k i (0, v)dA 



E ij (e,v)c iJk iE k i(e,v)de, 



(3.12) 



where v = ^nominal for short. It can be shown that the effect of infinitesimal changes of crack tip velocity (Sv) 
on the angular variation of the stress field is negligible and hence, the change of strain energy in the disk 1Z S 
due to the change of crack tip velocity can be simplified to read 11 



SAU e = i [k(v + Sv) 2 C Ue (H L ,a°°) ^ k(v) 2 C Uc (H L ,a°°)] ~ C Ue (H L , <r°°)k(v)k' (v)6v, 



where 



C Ub (Hl,<t c 



K 



10 



H L r s H 



E lJ (0,v)c lJkl ^ki(0,v)de>o. 



(3.13) 
(3.14) 



Fig. 3.2 schematically shows the behavior of the function g(v) = k(v)k'{v). Note that the above integral is 
positive and hence strain energy change is always negative; for all values of Cd/cu and v/cr strain energy is 
released as expected. 
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Figure 3.2: Schematic plot of the strain energy variation function, g(v) = k(v)k'(v) versus normalized crack tip speed v/cr. 



Now, we assume that the change of kinetic energy due to the released strain energy is dominant in a disk 
IZk with radius r k . Similar to the case of determining radius of dominance zone of strain energy release, the 
fracture quantum (ao) is the characteristic length in the problem. Assuming that r& = rk(ao,H, 

C-R-i ^nominal J j 

dimensional analysis tells us that: 

r k =a Q e(H,^^.), (3.15) 



CR 

where cr is Raylcigh wave speed, and is a dimcnsionless function. Now we have the following relation for the 
amount of kinetic energy change in the disk lZ k : 

AT= I -piiiUidA, (3.16) 



11 Note that A(*) denotes change in the quantity * due to crack growth (change in the crack length) while S(*) denotes change 
in the quantity * due to change in the crack speed. 
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where in is the velocity components of material particles. At the radial distance r from the crack tip, the particle 
velocity depends on the crack tip velocity and the strain at the particle position. The velocity of particles in the 
dominance zone of kinetic energy, TZk, is also a function of the angular position of particles and the ratio of the 
crack tip speed and a characteristic wave speed. If we assume that for the limiting roughness Hl the effect of 
the change of the nominal crack tip speed on the size of the dominance zone is negligible, the dominance zone 
of the kinetic energy change is fixed, and once again these conditions dictate that kinetic energy can be written 
as: 



AT = fe(t>)V 



1 



K 



-i 2 



dA, 



(3.17) 



where v is short for w nom inai- We know that the effect of infinitesimal changes of crack tip velocity (Sv) on the 
angular variation of the particle velocity field is negligible, and hence the change of kinetic energy due to the 
change of crack tip velocity can be written as 



1 



SAT = -[k(v + 6vj 2 (v + 6v) 2 CT(H L ,a 0O )~k(v) 2 v 2 C T (H L ,a 00 )] 
~ Ct(H l ,<t°°) [k{v)k'(v)v 2 + k(v) 2 v] Sv, 



(3.18) 



where 



K 



E 



2jt 



^i(e,v)^i(e,v)de > o. 



(3.19) 



Fig. 3.3(a) schematically shows the function f(v) = k(v)k'(v)v 2 + k(v) 2 v for different values of £ = Cd/cR. Note 
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Figure 3.3: a) Schematic plot of the kinetic energy variation function, f(v) = k(v)k' (v)v 2 + k(v) 2 v versus normalized crack tip 
speed v/cr. b) Terminal velocity versus £ = c^/cjj. 



that in plane strain Cd/cR G [1.635, oo], and for plane stress Cd/cR G [1.635,2.145] and hence the appropriate 
range of Cd/cu should be considered for each case. 

The schematic graph of the function k(v)k' (v)v 2 + k(y) 2 v in Fig. 3.3(a) shows an interesting phenomenon: 
After reaching some value of crack tip speed (vl) increasing v the change of kinetic energy becomes negative. 
Because the particle velocity inside the dominance zone of kinetic energy is proportional to the crack tip speed 
increasing the crack tip speed the kinetic energy change must be positive, and therefore the crack can not pass 
this limiting speed (i>l). In other words, when the crack tip speed increases the kinetic energy of the region 
around the crack tip must increase as well, i.e. kinetic energy must be an increasing function of the crack 
tip speed. In terms of the equivalent crack tip mass introduced by Chekunaev and Kaplan [2008, 2009], our 
argument is equivalent to saying that the equivalent mass must always be positive (note that their equivalent 
mass is independent of velocity), which is the case as was mentioned in §1. 

The predicted limiting speed for Glass, Homalite-100, PMMA, K5 (Glass), K6 (Glass), and SF6 (Glass) are 
as follows: 
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Glass 


(f 


= 0.220): 


I'L 


= 0.321c , 


VL 


= 0.320c , 


ii) 


Homalite-100 


(" 


= 0.310): 


VL 


= 0.311co, 




= 0.310c , 


111 ] 


PMMA 


(v 


= 0.350): 


u L 


— 305cn 




— 306rn 


iv) 


K5 (Glass) 




= 0.227): 


I'L 


= 0.320c , 


I'L 


= 0.319c , 


v) 


K6 (Glass) 




= 0.231): 


VL 


= 0.320c , 


I'L 


= 0.318c , 


vi) 


SF6 (Glass) 


[y 


= 0.248): 


T ! l 


= 0.318c , 


I'L 


= 0.317c , 



i ~ plane strain i - plane stress T , • r m i_i i 1. j_i j_i 

wncrc vl = v L and vl — . It is seen irom lable 1. tnat tncre is some scatter m tnc 



Material 


Author 


v l 1 Co 


UL I Cr 


Glass 


Schardin and Struth [1938] 


0.30 


0.52 


v = 0.22 


Edgerton and Bartow [1941] 


0.28 


0.47 




Bowden, et al. [1967] 


0.29 


0.51 




Anthony ct al. [1970] 


0.39 


0.66 


PMMA 


Dulancy and Brace [1960] 


0.36 


0.62 


v = 0.35 


Cotterell [1965] 


0.33 


0.58 




Paxson and Lucas [1973] 


0.36 


0.62 




Fincbcrg ct al. [1992] 




0.58-0.62 


Homalitc-100 


Bccbc [1966] 


0.19 


0.33 


v = 0.31 


Kobayashi and Mall [1978] 


0.22 


0.37 




Dally [1979] 


0.24 


0.38 




Ravi-Chandar and Knauss [1984a,b,c,d] 




0.45 




Hauch and Mardcr [1998] 




0.37 


K5 (Glass) 


Senf, et al. [1994] 


0.29-0.3 


0.5-0.52 


v = 0.227 








K6 (Glass) 


Scnf, et al. [1994] 


0.27-0.3 


0.47-0.51 


v = 0.231 








SF6 (Glass) 


Senf, et al. [1994] 


0.2-0.23 


0.34-0.4 


v = 0.248 









Table 1: Experimental values of limiting crack speeds for brittle amorphous materials. co,cr, and v are the longitudinal sound 
wave speed, Rayleigh wave speed, and Poisson's ratio, respectively. 

experimental data and this makes any comparison with experimental data difficult. However, we see that our 
estimates are close to the experimental data. It should be noted that what we have obtained for terminal velocity 
is an upper bound. 12 Interestingly, all the experimentally measured velocities (except PMMA) are smaller than 
our prediction of the terminal velocity. 

The relations between various wave speeds and Co are: a = \j (\+l)(\-.iv) c ® (pl ane strain dilatational wave 

speed), c p d = 2 ^ (plane stress dilatational wave speed) , c s = ^ (shear wave speed), cr = c s (l— 

(Rayleigh wave speed), where k 2 = ^(i-v) ^ or pl ane strain and k 2 = for plane stress. In Fig. 3.4 the 
normalized limiting velocities vl/cr and vl/cq are plotted for different values of Poisson's ratio and for both 
cases of plane stress and plane strain. The calculations show that depending on the Poisson's ratio the limiting 
velocities are in the range 0.276cq — 0.341co for plane strain and in the range 0.290cq — 0.341cq for plane stress. 



12 Note that any velocity larger than our calculated vl leads to a decreasing kinetic energy change that is not physical. However, 
this does not mean that a lower terminal velocity is not possible. In other words, what we calculate is an upper bound to terminal 
velocity. 
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Figure 3.4: Normalized limiting crack speed versus Poisson's ratio in plane stress and plane strain. 



4 Conclusions 

In this paper, we first obtained the asymptotic stress field around the tip of a dynamically propagating self- 
affine fractal crack. We then showed that there is always a lower bound to roughness exponent. We next 
looked at crack propagation and the asymptotic behaviors of kinetic and strain energy changes due to crack 
growth. We obtained an upper bound for terminal velocity by postulating that the kinetic energy change must 
be a monotonically increasing function of nominal crack tip speed. We predicted a material-dependent terminal 
velocity in the range [0.500cr, 0.557cr] and [0.539cr, 0.557cr] for plane stress and plane strain, respectively. 
We should emphasize that our asymptotic analysis only gives an estimate of terminal velocity. It was observed 
that for several amorphous brittle materials our predicted terminal velocities are in good agreement with the 
experimental data. In summary, our main results are: i) Fractal cracks tend to reach an approximately constant 
terminal roughness close to Hl = 0.8. ii) There is a terminal velocity of crack propagation lower than the 
Rayleigh wave speed, iii) The terminal velocity is a material-dependent fraction of the corresponding Rayleigh 
wave speed, iv) This material-dependent fraction only depends on the Poisson's ratio. 
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